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SUMMARY 



A theoretical investigation is made of a circular jot issuing 
from an orifice into an air stream flowing parallel to the jot axis 
and in the same direction as the jet. The analysis is aivioed into tv.o 
portions i 

1. The region immediately downstream of the orifice in 
which exists the potential cone of uniform velocity* 

2* The region downstream of the end of the potential oone 
in which the jet is fully developed. 

The solution for these two portions are fitted together at the section 
at which the potential cone disappears giving a complete solution for 
the whole jet* 

Numerical computations are made for several ratios of free stream 
velocity to jet velocity to determine the spread of the jet, velocity on 
the axis of the jet, and spread of the surface on -vdiioh the velocity is' 
one half the sum of the jet velocity and free stream velocity. 

Squire and Trouncer’s results are plotted as a comparison with those 
obtained in this analysis* 
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I. IKTROul'CTION 



ITie purpoGO of this analysis is to provide a conperison of the 
characteristics of a round turbulent jet disohorp;in{; into a noving 
air stream rrith the knovm characteristics of a jet discharging into 
air at rest. Both jet and free stream are assumed to have the aame 
density and the flow is assumed to bo incompressible. 

At the present time there exists no satisfactory theories for 
turbulent mixing. However, the integrated forms of the equations of 
motion, along with some dimensional reasoning, oen be used to give 

If 

some useful results. 

The solution of the round jet in a general stream has been 

already carried out by Squire and Trounoer, In their solution the shear 

( 2 ) 

stress is based on Prandtl's momentum transfer theory which is Ichown 
to be in error. It therefore seems worth while to investigate this 
problem using the integral equations and also making use of experimental 
data for obtaining the shear stress in the mixing regions. 

Use of the integral form of the equations was first proposed by 

(' 3 ) ( 4 ) 

von Rarman ) one of the earliest applications was by Polhausen , in 

( s) 

investigating laminar boundary layers in a pressure gradient, Sutton 
originally used the integrated form of the mechanical energy equation 
in laminar boundary layer investigations, uiepmann and Laufor*^iiave 
also applied this method to some simple two dimensional turbulent 
flows. 

In the integrated form of the equations of motion the only 
assumption neerssary is that for the velocity profile across the 
nixing region. This can be made quite accurately and for e closer 
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oompariaon to the work of Squire and Trouncer the se/re velocity profile 
that tftey use will be assumed in this solution. 

In transforminr, the nomenttua equation to tho inte.i?;ratoa form 
the term conteininr the shear drops out. It is therefore necessary to 
introduce the ujechanical enerj^' equation, 'fhis eq..ation, while not 
independent of the momentum equation, is necessary in order to introduce 
a term contain! njr the shear stress. 



J 



The folIoR'ing notation Trill to used. 




Kozzle Exit End of 

Potential Cone 



Fig. 1 

X Distance from noezle exit along jet axis 
r Distance from jet axis 

x^ Distance from nozzle exit to end of the potential cone 
a Radius of the potential cone at any point 

b Radius of the outer jet boundary dt any point 

b Diameter of nozzle exit 

u Velocity parallel to jet axis at any point 
Uj^ Jet Telocity at nozzle exit 
u^ Free stream velocity 
Uq Velocity on the jet axis at x —x^ 

Tjj Velocity of the stream normal to the jet axis 

Q Stream density 
T Shear stress 
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II, MALYSiS 

A, Region between the exit and the end oC the potentii^l cone. 
The three fur.daiaental equations are 



Continuity 



dCru) 

ax 



+ 



a(ru) 

dr 



= o 



Momentum 



QU 



dm 

ax 



-j- py fr?^) 

^ dr . r dr 



Mechanical energy 



PU.^ +puif ^ - -!=L -2-rrTr') 
dK ar r ar'^' 

These three equations are integrated between the 
boundaries of the mixing region and are then 



C ontinuity 
f b 

urdr +u-^,b + ^ 



d 

dX 



wa 



Momentum 



~Ufcb ^=0 



A 

dx 



va 



^uADtu:aft-utb^ 



uVdr +UuVv,b 



-0 



Mechanical energy 
r>b 

U . M. . 2 

a 



I die 

Idx 









a 



A function is assumed for the velocity profile 
across the jet such that 

U = dfe+ - COS.TT^] 



or 






whore 



dX. 



^ b-a 

r = b-(b-ci)'^ 

dr--(p-a) 6 i 

Substituting into the continuity equation 
ro 



Uk+ 






%)][b-(b-a)d:^ ][-|)-a)d:^] 



+Vtb+u«,a ^ -Ufcb =0 



Expanding 

_d 
dx 



Ujab-b")|d;? +Ubfa-br|^d5? ■^^^Yab-Wp(^)d^ 



+ ^(a-bf 




+7/^b + Ua.0 Ufcb —O' 



Integrating 



A 

dx 



-u^(ab-b'*) - Y(®"bf + (ab-b^)!, 



+ ^^Va-bnJ +Ubb-(-U«a4| -dbb^ =0 



where 



I, = l\k 7 i)dyi 



6 



Dif I’erentiatins: aud collectini? terms 

[AQ +Bb] ^ + [B-a + C b] 41 -0 



(I) 



where 



A — (Uq.-U.(p){I 2+ 0 

B = (u^-Ub)(% - I2) 

C =(Ua-U,)(l^-I,) 

Substituting into the momentuni equation 



+ \X^\yJo +iXla ^ -u,^b 41 - 



0 



Expending and integrating 

^db-b J +■ t^6xb~b^) I, 

+ Ufcir^b -vala 



da. _ 

K 



a.^b 44 - = 0 - 



where n 

I. . I 

Lifferentiating and collectinf terras 

fOa + Eb] ^ -t-[Ea+Fb]^ i- u^Ut b =0 ( 2 ) 
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where 

D— ^ + 2. T2. 

E ~ C'^ ~ +Cl,~2ij)Ub(Ui-Ub) 

F = (^- ^(^^-Wk)" + (2i,-2I)a^(u,.-uO 

The shearing; stress now appears in the mechanical 
energy equation and an assumption will be made for the 
shear profile* The shear stress is a variable and could 
be determined from the velocity profile . However, a 
shear profile derived from an already assumed velocity 
profile would have u greater degree of error due to 
small errors in the assumed velocity profile. For 
simplicity an average shear profile will be assumed 
for the whole region. 

From experiments with a jet discharging into still 
air shear profiles have been found as shown in Fig, 2, 

At the jet exit the shear profile is symmetrical. Past 
the end of the potential cone the maximum shear occurs 
at about '^=.667 to ,75 and the shear profile has a steep 
slope on the high velocity side. An average shear profile 
for the potential cone region is shown in Fig. 2 and is 
given by 

r =-e 

which is assumed to be of the form 

r = ^(Ua.-UtfgC^) 



where 



g ('>1) -- 



a'v' 



(u^-u) 

In particular it will be assuraed here that the 
'unction can be represented by 

g ^ cos"^ Vi C^~ 



Substituting into the mechanical energy eouation 

.... . -t3r 

UuH- 



z 3 T^ 









I 






_ _ _L 

~ Q 







Expanding and integrating 

A‘ 



-Ub ab(a^-<u^)I,d- 



-h 






Cu.-u,)%][ab-b^ 



-t- 



+ ^ Ub - u,y-l I J [a- b] S 

+ Ub^TTbb + - <^k.b 

= TT ^Uo.-Ubf (IT--Ig)b + TT (U,^-Ub)^l 8 a 



where 



Is= 

Is = 

I. = 

Ia= ff'C^)^Cp 
Differentiating end collecting terras 



Ga 



Hb]^ 



+ 



Ha -f-Jb 



db 

■ar 



-bU^lTub = Ka+ Lb 




where 



G = a/-o/ +3^^Cu^-u^,)I^ -t-^uj.u^-uSl4 

+ ^ Is 

M= + (/(b(W.-v- J^3~l J 4 ) 

^ (.^K “ C 8 I 5 ~ 

J= U,"(^^4,-WiX31,-3I,) +Ub(i/^-6/jYil4-|l) 

-t- Cu.4- - ?Is) 

K= ■n' (UoUb)® Ig 
L= TT (u^-uj(l,-lg) 



Eliminating in equation (1) and equation (3) 

Ga + Hb] ^+[Hq + Jb] ^ [Aa +Bb] ^ 
-Ub[bq+Cb]-^ = Ka + Lb 
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\ 



+|j“UteCj'^ — Kq + Lb '■ 
d r 

^ [ cxaSz^ab + 'zrb^j = Ka + Lb 

■where 



oc - &~~ ^ _ H ~ Uu B 

a ? z 5 

Using the chain rule of differentiating, 
the above equation becomes 



^ = J~^b c 

a 

cJX.) _ di ) , ^ 
d cl b dx 



= l\a + L D 



[Zo,a^ -^2'3'b]^ 

(Ka + Lb)-^ = Czota + 2(3 bj) ^ +Zpo. + 2'3’b 

Eliminating in equations (1) and (2) 



Ba-bEb 



dcx 



-fjEa-Hpbj 



u. 



dib 

j dX 

j^B(X "hC bj 



u. 



/Ka-v-Bbj 



d (\ 
Ix" 



^ -o 



+• [E-Ub&][a 






dx 






(4) 



db _ 



dX 



- 0 



Integrating 



Sa^+ 2 eab +• 4 ^ i-Z =o 

where 

8-=5^i^) e=£r^^ - P-u^,C 




2^ — constant of integration 
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Dlf forontiatinis; with respect to b 



2Sa§|+2€a + 2eb^ + 24>b=0 



dx _ _ eg -V 4>b fe) 

d b 6 b ^ Sol 

Substituting equation (6) into equation (4) 



jka+Lb ^ 



■V2 



P ^ X - ftrr ?' 



2 (3a +2rb 



1. ^ -<^6)a'^^-(V^ ~ g<t')gb b 

2d(o Gb“bScL 



1 djL _ \ Az b As b 

2 db ' So. + €b Kcx+Lb 

whore 

X , M. z..5g_ 

5K 

<S(eK-5L) 

~ ~ boc^ ~t- ~ <^(^) 

K(Lb“K€) 



From equation (5) 



5 a .+€ b = +- b ^(€^-0 6 )' 
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■JL dji r A + ^ ^ A3 b 

^ Ka+Lb 



Inte grating 



X _v L +bYe"-4>S) . 



bdb 



K<id" ^ b 



This dxprcsslon oah b6 in^6gr&t6d and gives 

|^=x,b- \:€HEEEi^ 

^ 4>6-€^ 



'“9 ^ 



c. 1^^ b^(C;.-M' -c .- TCc^'nXc.Vc,'? !' <• c, 

? + cA(et+iyc,’--<:^b' - c, 



+Y 



«^er« 



Co = 






LS-K6 



P 2 _ _ K'^gX 

■ ClS- Kt)"^ 



C. = 






(L£-^€.)^ 

Y * constant of integration 
Y can be found by the boundary condition that x =0 
whore b = d/g. Obtaining the parameters x and b in 
dimensionlesfl form the result becomes 



X /b A eA. 



+ ZX-t-/- ('t>&-6^J' 



4>S -e^ 



+ A 3^ o rb 

Cz+ 






c. 



a-fciTi 



log 







1'.*' ’ a 


^('c^+D 





XiC„ I _yfz^ ' 4 _ 

c;?rL' ^ " 






c, I &c^+i' - 

c.] - ^'J 



Tho constant of integration ^ is obtained from the 
boundary condition that a = b = <^2 x = 0 so 

J’i. 

Z. - — -g- U«, (u.a.-u.b) 



(7) 



The only rtnaininr constant to be detemined is 

which appears in the shear term* Prom experimental data 

^ f6) 

obtained by Liepmann and Laufer' it is found that near the 

jet exit ^ ~ 0,0072, From experimental data 

(7) 

obtained by Corrsin it is found that in the fully developed 
region ~ 0*022, Both of these results are for 

the case of a jet discharging into still air or ^ =0 
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Differentiating our assiuaed function it is found that 
occurs at the point -where ^ = 0,692 and ~ ^ 

The constant 'h must be given some average value such that 
the average shear profile -will be bet-ween 0.0072 
and 0.224. 

( 8 ) 

From experimental data it is known that the end of the 

potential cone is at a point about 9 radii downstream of the 

Uu 

exit of the jet for the limiting case where u^ = 0. In order 

X. Ci 

for “ ^ =■ 9 vhBxi O and a = 0 

and "fe = 0.07B88 



BOff 



B. Region down«troeni of th® end of tho potential cone. 

At tho point dovmetream of the jet exit where x 
tho potential cone diaappears and the rariable, a, no 
longer appears, DownstreaLin of this point the velooity on 
tho jot axis, Mq , will become a variable decreasing. with x. 
The assumption for the velocity profile in the region of 
the potential oone was 




At -= X - a - 0 and u ^ = u 




This roduooo to 




or 





whore 



The continuity aquation in this rep;ion than bacoraes 



A 



r 






4 u;b-abb^ = o 



This can be raducad to 
Trtiara 

I 9 = 

V = Uc-u^ 



. Substituting in the momantunj equation 




Ub(u^-Ub)b"'j^^-f(yU.)C^ 




+ 



ju.-f (ju-)d^ b ^ -0 



« 

^ V?hich raducas to 



[4Isu,V+I,„V^^ 




2u^,irb = o 



i*hara 



^10 



17 



In the mechanical energy equation the shearing etreaa is 
given by 

The function <^(jx) must be such that tho shear profile will 
resemble that illustrated in Fig, 2, Tnis will be 

= I J, ^ -f djx 

Z = constant 



Substituting into the mechanical energy equation 

^ [Ub b'" I + I Ub\u, - Uj,)b^J' 



which reduces to 



Ufe V ■+ +4l„V 

+ [|I,u^b +1 1,„UbbV +i I„bV^] 

whore 

1 1. ={/xlV)d>J^ 




(loj 



Eliminating in equation (8) and equation (9) 

[4Uu,V + Lzi,u, + l„v]b 
+ 2a,[-l3V^'-iI,b4^] -0 

db _ IgUb + I.oV olV 

b 2l9Ufe + I,oV V 

Integrating 

log b =-^lo 3 [^IgU^-tl,^v]-ilogV 

cl _ b -bX/oV ^ d V 

~ ~ V 2XgU^+l,<,V ’ V 
d b _ _ C,Cl &Ub **• IiqV) dV 

d7 ' v'^"(zi9Ub+i,oVf" 



1 



Elimineting in equations (0) and (10) and aubstitutlnf; 
for ^ above we rot a single equation which is 



Islio^b -V- I glu ^ 



dv ^ 



2_Tji. \j 

C, 



This my be integrated to give 



8 In. le , j y 

^ j -r ^ 

-L \\ 



-ii- 

31.. 




V;> 






The constant of integration C jnay be found by the 
boundary condition that at x=x^, 




The constant of integration C, may be evaluated from 
these same boundary conditions and from equation (5) when 



a= 0 
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Subatituting those values of C, and C 2 ^ the final 



solution is obtained. 

r 



3 ? 2 . 



tet- ^ ^i.oO- 









'/z 



J 



X 



If- 


- 

C 


Vi. 

Ir. 


< 

- 1 - 


~1 T T - T 

^ * lO -L U / 0 




K 


r r J 


61 ^ 1 ,, J 






( 13 ) 



where 



% 



Uk +V 






V 



Equation (11) becozsas 

b' - 

1) 



-1 


[ fs - a '- ft ] 


|[ l , o ^ k +(^ l 9 - I.oi 


1 









(l4) 



Xj. may be found by solving for b when a=0 in 
equation (6), Using this value of b find the corresponding 
value of X in equation (7) which will be x^. , 

Equation (13) gives x os a function of (u^-Uj^, 
Equation (14) gives b as a function of (u -u,). It does 

Co 



not seem necessary to get an explicit expression for x 
and b. By assuming ralues of (u^-u^) correspondlnrv values 
of X and b are found and b is plotted against x in 
Fig. 5. 

I\OY/ the constant Z appearing in the shear stress 
expression must be evaluated. In the fully developed region 
downstream of the end of the potential cone ^ 0.022 . 

From the assuaied function of ^ (/a.) it is found that 
= 0.4382 or Z =0,05116. 



IIIo DISCUSSION 



'i2 



The roiults plotted for the region containing the potential 
cone are the results obtained by Squire and Trouncer* In using the 
derived equation for x and b (equation 7) in the potential oone 
region the final results obtained did not appear physically reason- 
able, particularly at higher velocity ratios of free stream to jet. 
Due to the enormous amount of computation necessary to recheck these 
results Squire and Trouncer *s results were used in this region. 

The very close ccrinparison of the two analyses in the fully 

I 

developed region would seem to indicate that the integrated equation 

t 

solution that has been used here should give fair agreement in the 
potential cone region. 

The basic analysis of the potential oone region is known to be 
correct through equation (5), This equation gives a relation between 
the spread of the jet and the width of the potential cone. The 
remainder of the analysis within this same region or the computations 
may be in error. However, the method of solution is outlined for any 
future investigation of this problem. 
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AP?I7fDIX I 



Evaluation of dimensionless integrals appearing in the 
equations of motion. 

The dimensionless integrals appearing in the equations 
of motion all may be expanded into integrals of forrj found 
in standard integral tables. Following are the results of 
these integrations, 

I, = [ =|[l- costt)^]oIt^ = -1,00 

=1 [i-co5Tr^J^d)2 =-a^02 
Ijf ] -f (gjd>^ =JJ)-co5Tr:^]^ol>^ = - i.soo 

= Jjl- cos TT)^]\d)^= -1.155 

^s'l ~ I (^1 'Cos IT- d)^ = -Z.SOO 

^6 "i I = -2.014 

gC'>()T(>^)<i>^=J'fe)^^Cos’Tr(5j-i)si>ni',^d:^ =-. 0687sx(o'^-fe 



I.= 





*'eo5^iT('?^-4:)S)nfr>jc/)^ _-,043osx^o ^ 







p + Co6TT^]yud/A = 0.298 



1,0 “ j' [\ f CostTjjII''m-‘^J^ - 0.346 





! +COS rr/M /ac/ju. = 0.4-85 



/U,g(yLg-fpA.)cpA. 



fl'Isin iT/Ajj + tosT\jx ^ 4- -L (cos Tr/A,i-TT)A.sin rrju. - 



= -0.03982 



APPENDIX II 



Integrations necessary for solution of the equations. 



A, Integral appearing in the equation for the region 
containing the potential cone. 



1 = 



b cfb 



bCoL'l" L b 

Substituting for a froa etuation (5) 



1 = 



b db 



Keb i-K\/4€^b"-4S(»b"-£)' + 1 u 
8 



which retiuoes to 



n 



I = Co 



b db 

b Wc,’'- c^b^' 



where ^ 

r 2 

° LS-K6 

K^£S 

' CL&-K€)^ 

r - -€"') 

CLfe-Kfr)"*- 

Rationnlieing the denominator 

f (b^- j d b 

b'^(c^t-i) - c,^ 



I = Co 



c, 



b^'clb 






- Co 



Vc,^- b d b 
b^(c^+i; -e 



d b 



b^(c^-t-i) -c 



-1. 



b^d b 



[_b^/ C2.+T' - C,][bV?2.^-^ 



n 



C, 



C . 



+ 



2Cc*.+ l) 



Z(C^i) 



bVQ+? -C, bVc,+-l'4-C, 



db 






C, 



C, 



bN^cJ+T'C, b^/c7T^ i-c, J^(C^+/) 



db 






Z(q 









" 42 I 5 E; bdb 



bYc^->-i)-c’ 

Le't c,^-Ci.b = CiRo 

~Zc^bd>b= ZC2_RodR^o 

bdb = - RodRo 






c.'-c.R,'- 






■/T, Ro(-~Rocn?0 

■ J Cc.+i)(-|^ --r:) -C,^ 
■VK RqMr, 



_ ^ 

2(C:,-t-0 



2 + 



-Jcl 






RoVc^+i' - 



de, 



VcT 

2 Cc,+ i) 



l_ 

2(C^-H) 



21^ ^ c, I03 (Ro^g^M' - c, icx^(Ro^ic^+^-^-^ 

Cc^i-|) ^l C,(c^+|) . 



2V2:^ + 7 ^ log 



B. Integral appearing in the equation for the region 
downstream of the end of the potential cone. 




dV 






dt 

(l-t-)- 





I.O 




± 



4 



eft 





Lls±. 

3 I, 




c VI^ iJIBMi jJIM ^ r»v 





J 

•V 

< 



c 











a 

<v 



n. 



> 

0 

1 





rr 

$ 

i 

K 

n 

O 



D 

< 



O 

•n 

< 






I 
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